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ABSTRACT 

Context. The investigation of the pulsation properties of pre-main-sequence intermediate-mass stars is a promising tool to evaluate 
the intrinsic properties of these stars and to constrain current evolutionary models. Many new candidates of this class have been 
discovered during the last decade and very accurate data are expected from space observations obtained for example with the CoRoT 
satellite. 

Aims. In this context we aim at developing a theoretical approach for the interpretation of observed frequencies, both from the already 
available ground-based observations and from the future more accurate and extensive CoRoT results. 

Methods. To this purpose we have started a project devoted to the computations of fine and extensive grids of asteroseismic models 
of intermediate mass pre-main-sequence stars. The obtained frequencies are used to derive an analytical relation between the large 
frequency separation and the stellar luminosity and effective temperature and to develop a tool to compare theory and observations in 
the echelle diagram. 

Results. The predictive capabilities of the proposed method are verified through the application to two test stars. As a second step, 
we apply the procedure to two true observations from multisite campaigns and we are able to constrain their stellar parameters, in 
particular the mass, in spite of the small number of frequencies. 

Conclusions. We expect that with a significantly higher number of frequencies both the stellar mass and age could be constrained 
and, at the same time, the physics of the models could be tested. 

Key words, stars: variables: 6 Scuti - stars: oscillations - stars: pre-main sequence - stars: fundamental parameters - 

1. Introduction out that the interior structure of PMS stars entering the in- 
stability strip differs significantly from that of more evolved 
Pre-mam-sequence (PMS) 6 Scuti are intermediate-mass main-sequence stars (with the same mass and temperature), even 
(1.5Mo < M < 4M0) stars that cross the pulsation instabil- ^^^^^^ structure of the env elope is sim ilar. This property 
ity strip of more evolved classical pulsators during their evo- subsequendy confirmed bv ISuran et all JMH) who made 
lution toward the main sequence (MS). The first evidenceof ^ ^ comparative study of the seismology of a 1.8 Mo PMS and 
Scuti pulsation in PMS stars was provided by iBregeij (11972D, po^^.^s stai- and found that the unstable frequency range is 
who discovered two candidates in the young cluster NGC roughly the same for PMS and post-MS stars, but that some non- 
2264. More than 20 years elapsed since this discovery before ^^^-^^ ^^^^^ ^^^^ ^^^^-^^^ ^j^^ ^^^p -^^^^^^ 
other authors confirmed the occurrence of this kind of pulsa- j^e star. In particular, it is possible to discriminate between the 
tionobserving the Herbig Ae stars HR5299 d Kurtz & Marand , pMS and post-MS stage using differences in the oscillation fre- 
1995|) and HD104237 jPonatietaU [19971). In the same years q^^^^y distribution in the low frequency range (g modes, see 



Marconi & Pallal (I1998D obtained the first theoretical instability also lTempleton & Basull2003h 



strip for PMS intermediate-mass pulsators by computing an ex- 

tensive and detailed set of nonlinear convective models for the The theoretical analysis by iMarconi & Pallal (Il998h stim- 
first three radia l mode s, along the PM S evolutionary tracks by ulated an increasing observational interest for this class of 
iPalla & Stahled (Il993h . More recentlv iGrigahcene et alj (l2006l) pulsators and the cuiTent number of known or suspected can- 
have produced a theoretical instabily strip for PMS stars for the didates amounts to about 30 stars. An updated list is available 
first seven radial modes. Marconi & Palla (.1998.) also pointed at http : //ams . astro . univie . ac . at/p rns_co r ot . php , and there 

are re cent reviews by Zwintz et alj ( |2004|) . iMarconi & Pallal 

Send offprint requests to: A. Ruoppo (|2004|) and .Ripepi et al.. (120061 . However, only a few stars have 
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been studied in detail, so that the overall properties of this 
class of variables are still poorly determined. In this context 
our group has started a systematic monitoring program of 
Herbig Ae stars with spectral types from A to F2-3 with the 
following aims: 1) to identify the largest number of pulsating 
objects in order to observationally determine the boundaries of 
the instability strip for PMS S Scuti pulsation; 2) to study in 
detail through multisite campaigns selected objects sh owing 
multiperiodicity (see [M arconi et al. 2001; Ripepi et al.l 120021 
I2003L 120061: iPinheiro et al...2003; .Bernabei et al., .20041) which 
are the best potential candidates for asteroseismology. 

In order to establish a theoretical approach for the in- 
terpretation of these data, as well as to prepare the basis 
for the interpretation of future more accurate satellite ob- 
serva tions expected from CoRoT (iBaglin & the CoRoT TeanJ , 
I2003h . in the context of the CoRoT/ESTA collaboration 
yittp : //www . astro . up . pt/corot/ ) we have started a theoreti- 
cal project devoted to the computation of asteroseismic mod- 
els of intermediate mass PMS stars covering a wide range 
of stellar masses. In this paper we present our first results 
based on an extensive grid of ev olution a ry PM S models com- 
puted wi th the CESAM code jMorell 1 19971) . with updated 
physics jMarques et al.L |2004|) relevant for PMS modelling. 
The Aarhus linear adiabatic non-radial pulsation code ADIPLS 
(http: //astro. phys.au.dk/~jcd/adipack.n/) has been ap- 
plied to the models of the grid in order to build a reference base 
to reproduce the observed pulsation frequencies. 

The paper is organized as follows: in Section 2 we present 
the evolutionary PMS models, we introduce the adopted non- 
radial pulsation code and report on the predicted pulsation fre- 
quencies; in Section 3 the computed large separation and its de- 
pendence on luminosity and effective temperature are discussed; 
in Section 4 we present a method to reproduce the observed fre- 
quencies and a theoretical test of the predictive capabilities of 
the method; in Section 5 we compare the theoretical results with 
the observations from multisite campaigns of the stars V35 1 Ori 
and IP Per. Finally the Conclusions close the paper. 

2. Reference grid of models with pulsation 
frequencies 

For exploring the pulsation properties of intermediate-mass PMS 
stars, we decided to build an extensive grid of models and fre- 
quencies to provide a reference on the seismic characteristics of 
PMS stars across the HR diagram. 

2.1. PMS evolution models 

The stellar models have been computed with the CESAM code 
dMorel, 1997) that is optimized for asteroseismology. The pre- 
cise calculation of oscillation frequencies (p and g modes) re- 
quires numerically precise and adequately dense meshes with 
the quantities describing the internal structure and entering the 
oscill ation equations (e.g lChristensen-Dalsgaard & Berthomieul 
1199 lb . The physical and numerical pro perties o f thes e models 
are discussed in detail in the paper by Morel (19 97|) with an 
update d version of the physics being used as in .Marques et al.l 
(l2004l) . 

From the CESAM PMS evolutionary tracks (see Fig. |3) we 
selected 56 models for stellar masses varying from 1 .6 to 4.0 Mq 
with a step of 0.2 Mq. For each selected mass, we consider from 
two to four different effective temperatures in the range 6 000 to 



10 000 K, along the corresponding evolutionary track. We report 
in Table [T] the physical properties of the selected PMS models, 
which cover the empirical HR diagram location of known pul- 
sating PMS intermediate-mass stars. 



Table 1. Physical properties of the selected PMS models. 
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Fig. 1. Predicted frequencies as a function of the radial order n 
for models at 3.2 Mq (bottom panel: mod33 left, mod36 right), 
and 3.6 Mq (upper panel: mod43 left, otoc/46 right). Symbols are 
connected by a full line in the case of l-Q, by a dashed line for 
l-\ and by a dotted line for 1-2. The effective temperatures of 
the models are indicated for each mass. 



2.2. Calculation of the oscillation frequencies 

In order to explore the pulsation frequencies that characterize the 
selected PMS models, associated to both radial and non-radial 
modes, we applied the ADIPLS code. The pulsation frequencies 
are obtained as solutions of a set of linear equations for small 
perturbations (for details see lUnnoT t al. (1979)). 

Tlie physical and numerical assumptio ns used for the 
compu tation are discussed in detail by C hristensen-DalsgaardI 
(|1982|): [C hristensen-Dalsgaard & Be rthomieul (Il99l[) : 



IChristensen-balsgaard & Mullani(il994i) . 

In this work we concentrate mainly on p modes, i.e. acoustic 
modes, with a spherical harmonic degree / = 0, 1,2 and radial 
order n = 0, 1, 20. We only consider very low degree modes 
because these are th ought to be more likely detectable (e.g. 
[Oziembowski 119771 : thristensen-Dalsgaard & Gough 198!^. 
From a first inspection of the results we note that the predicted 
range of p-mode frequencies moves toward higher values as the 
stellar mass decreases, at fixed effective temperature, and as the 
effective temperature increases, at fixed stellar mass (see Fig.[T]). 
This behaviour confirms the well known relation between the 
pulsation frequency and the stellar mean density. 



Fig. 2. Large frequency separation as a function of frequency for 
the model with M = 2.0Mq and T^^ = 8184 K, for 1 < n < 18. 
Symbols are connected by a full line in the case of 1=0, by a 
dashed line for 1=1 and by a dotted hne for 1=2. 

rR I 

between modes of the same order, where Av ^ ( dr/c) oc 

{mIB?^^^^, with c the sound speed in the star. 

In Fig. |2]we show for model modS the behaviour of the large 
separation as defined in Eq. ([T]i, for / = 0,1,2, versus the fre- 
quency. 

On this basis a relation can be found between the large fre- 
quency separation, the luminosity and the effective temperature 
(the stellar mass is related to the luminosity for each effective 
temperature). For this purpose, we have computed a large sepa- 
ration Av in the asymptotic region as the mean of all values for 
5 < n < 20. 

A linear fit for all models in our grid gives the following 
relation: 

logL = ( - 12.81 + 0.02) + (4.35 + O.OS) log - 

-(1.681 +0.015) log Av, (2) 

plotted in the HR diagram in Fig. |3] for selected values of Av, 
ranging from 20 /zHz to 70 yuHz. This relation is useful to ob- 
tain an expected range of the large separation of an observed 
pulsating star once its position in the HR diagram has been esti- 
mated. In the next section we use this approach as a first step to 
characterise a PMS 6 Scuti pulsator with the adopted non-radial 
pulsation models and the observed frequencies. 



3. The large frequency separation 

The pattern of regular separations between frequencies of modes 
of different degree and order is a characteristic feature expected 
in /?-mode stellar oscillations. The frequency separations depend 
on the internal structure of the star so by measuring them we gain 
valuable diagnostics on the stellar interior 

In particular, in the asymptotic regime, it is convenient to 
use the large frequency separation defined as 

Av„,/ = V,,,/ - v„_i,/ - Av , (1) 



4. How to reproduce the observed frequencies 

In this section we present a methodology to compare the ob- 
served pulsation frequencies with theoretical ones. This compar- 
ison allows us to obtain information on the stellar structure and 
physical parameters, such as the mass, the radius, the effective 
temperature and the age. For simulating typical available data 
for observed stars, we "theoretically built" two PMS test stars. 
The first one, St arl, corresponds t o a model computed using the 
STAROX code dRoxburghl l2005h as determined for Task 1 of 
CoRoT/ESTA, while the second one, Star2, was obtained with 
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4.2 4 3.8 3.6 



Fig. 3. HR diagram showing PMS evolutionary tracks computed 
with the CESAM code for the labelled stellar masses (in Mq) 
The linear fit given in eq. ^ is represented by dashed lines, 
defined for a constant value of Av (as labelled, in yuHz), from 
20 fiHz to 70 jt/Hz. 



the CESAM code (Morel', '1997') after integration of the birthUne 
from iPalla & Stabler (1993) as the initial configuration for the 
evolution calculated by CESAM (IMarquesl 1200 6). The parame- 
ters characterising these models are given in Table |2] 

We adopted two different evolutionary codes for producing 
the two test objects (also changing the initial condition in one of 
these) because we need to ensure that the use of the same code 
for building the test cases and the grid of models used here to es- 
timate the stellar parameters is not the dominant factor determin- 
ing the success of our inversion procedure. The physical proper- 
ties of these two test stars are presented in Table |2] whereas the 
corresponding theoretical frequencies, for both radial and non- 
radial p and g pulsatio n modes, were computed using the POSC 
code (lMonteiroill996l) - again a different code from the one used 
to calculate the frequencies of the reference grid. The values of 
the frequencies, the spherical harmonic degree / and the radial 
order n, are reported in Table |3] 

Table 2. Physical properties (exact values) of the two test stars. 
These values were only known to one member of the team. 
The quantities provided to the other team members are listed 
as 'Observed' values. 



Mass 


Radius 


Luminosity 


log(reff) 


Age 


M/Mo 


R/Re 


log(L/Lo) 


K 


My 


Exact values: 










Starl 2.0 


1.862 


I. mi 


3.9251 


0.00 


Star2 2.4 


3.482 


1.5915 


3.8887 


3.72 


'Observed' values: 










Starl 




1.19±0.82 


3.92±0.06 




Star2 




1.58+0.07 


3.88±0.05 





Table 3. Frequencies, spherical harmonic degree, and the radial 
order as computed using the POSC code for the two test stars. In 
the last column the 'observed' frequencies provided as observ- 
ables for the two test stars are also given. 
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166 ± 3.5 
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188 ±3.5 
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196.8293 
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196 ± 3.5 
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224.3726 
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224 ± 3.5 
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7 
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265.4323 
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266 ± 3.5 
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268.3870 
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268 ± 3.5 


fu 
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293 ± 3.5 


/l5 


296.9251 
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297 ± 3.5 



We note that most of these frequencies have been verified to 
correspond to modes that are actually unstable, through the use 
of nonadiabatic computations. We also notice that only one au- 
thor knew the real stellar parameters and the mode parameters of 
the frequencies, while the others tried to find the best solution re- 
producing the observed stellar parameters and frequencies, just 
as if they were true observed data. In particular, the author that 
knew the properties of the two test stars has simulated real data 
by shifting the frequencies, the luminosity and the effective tem- 
perature relatively to the exact value and quoting an error not 
always compatible with the original (exact) value. The frequen- 
cies for g modes have also been added in one case. All values, 
exact and simulated observations, are reported in Tables |2]-[3] 

In both test cases the number of frequencies is high (10 or 
more) and with a few of those already in the asymptotic regime. 
Present data on PMS 6 Scuti stars are not so rich (only a few 
frequencies) but forthcoming observations from CoRoT are ex- 
pected to provided several tens of frequencies per star To illus- 
trate the present situation we consider later an application to two 
real cases observed in multisite campaigns. 

Based on the 'observed' data we have tried to establish if 
we could determine the basic stellar parameters. In doing so an 
approach has been developed in order to best use the available 
data when there is no mode identification and poorly constrained 
stellar parameters (as is often the case for single PMS stars). To 
summarize our approach the steps we follow are: 
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- combine the grid of models with the known luminosity and 
effective temperature to determine a range in mass and ex- 
pected large frequency separation; 

- estimate the large frequency separation from the frequency 
data; 

- reduce the mass range by using the observed large frequency 
separation; 

- use a few models within the possible range of mass and large 
separation to reproduce the frequencies in the echelle dia- 
gram; 

- based on the best representation in the echelle diagram able 
to reproduce the observations we provide the preferred stel- 
lar parameters and a tentative mode identification. 

This sequence is motivated by the need to use sequentially the 
most robust information contained in the frequencies. This is 
done by reducing as much as possible the effects on the anal- 
ysis of the unknown physics that is not included in the reference 
grid of models and frequencies. 

The proposed sequence of steps is executed for our two test 
stars as described in the next sections. 

4. 1 . The HR diagram 

On the basis of the reported range of luminosity and effective 
temperature we show in the Figs |4}{5] the predicted HR diagram 
position of Starl and Star2 respectively. In the same plots we 
also display the set of PMS evolutionary tracks from CESAM 
for stellar masses from 1.6 to 4.0 solar masses, together with 
the linear fit obtained in the previous section for constant large 
separations (ranging from 20 to 70 yuHz). These plots allow us 
to define the range in mass consistent with each star and deter- 
mine the possible range of values for the large separation. The 
obtained mass and large separation ranges are: 

Starl: 1.6 < M/Mq < 3.4 and 20 fiHz< Av < 80 ^Mz, 
Starl: 2.2 < M/Mq < 2.8 and 20 a/Hz< Av < 50 juHz. 

These are yet large ranges of possible values which correspond 
to the uncertainty on stellar mass when only the classical observ- 
ables (luminosity and effective temperature) are available. 

4.2. The large frequency separation 

In a real case the large separation would be ideally extracted di- 
rectly from the power spectrum of the light curve based on the 
expected range reported above. The development of an extrac- 
tion method for this quantity from the power spectrum, based on 
the methodology developed for solar-type stars but adapted to 
6 Scuti PMS stars, will be described elsewhere. 

For the test cases we have here, there is not a power spectrum 
so an alternative analysis has to be used in order to estimate the 
large separation from the actual frequencies. Based on the reg- 
ularity of the spacing between frequencies we obtain an indica- 
tion of the large separations as being about Avi ~ 65 //Hz and 
Av2 ~ 27 /iHz for Starl and Star2, respectively. This approach 
does not provide an uncertainty but in a real case this 'error' 
could be as high as cr(Av) ~ 5 juHz. 

4.3. The echeile diagram 

As it is well known, the echelle diagram is convenient to illus- 
trate in detail the properties of the frequency spectrum. In this 




4-2 4 3.8 3.6 



Fig. 4. HR diagram position of Starl (dotted square) together 
with CESAM PMS evolutionary tracks (full lines, labelled with 
masses in Mq) and the linear fit (dashed lines, labelled with 
Av in fjHz) of constant large frequency separation obtained in 
Section[3] 




Fig. 5. The same as Fig. |4] but for Starl. 

diagram the frequencies are reduced modulo Av by expressing 
them as 

Vn,i = vo H- A: Av -H V , (3) 

where Av is the large separation, vq is a suitably chosen refer- 
ence, and k is an integer such that v is between and Av. The 
diagram is produced by plotting v on the abscissa and vq + khv 
on the ordinate. If the asymptotic relation between frequencies 
and large separation were precisely satisfied, we would obtain 
data points arranged on a set of vertical lines corresponding to 
the different values of /. 

By using the best values for the large separation, determined 
before, we select all the models for which the intersection of the 
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Fig. 6. Echelle diagram for Starl, using the theoretical frequen- 
cies computed with the Aarhus code. The parameters used to 
produce this plot are shown in the plot. The observational un- 
certainties would correspond approximately to the size of the 
symbols representing the observed frequencies fi. 
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n degrees (see Fig.|7]i. The only exception is /i that corresponds 
to a g mode and this is the reason why the echelle diagram has 
been unable - correctly - to adjust this particular value, even if 
our best-fitting model has a g mode that fits fi within 0.005 fiHz- 
These are both ideal test cases, but they demonstrate that if 
enough frequencies are known with sufficiently high precision, 
the oscillation spectra can be interpreted using as reference a de- 
tailed grid of PMS models and their frequencies as proposed in 
this work. The key assumption is that the measured values of 
the frequencies do not deviate strongly from the predicted val- 
ues and relative spacing. Any effect that changes significantly 
the linear adiabatic radial and non-radial frequencies used here 
as the reference may invalidate this assumption. These can in- 
clude fast rotation and strong magnetic effects on the frequen- 
cies. Before a definite characterisation of the observed modes is 
reached, these effects need to be taken into account on a case- 
by-case basis. 

The heavy-element abundance Z (or metallicity) is also an 
important parameter of the grid (here Z - 0.02 has been used) 
that has an impact on the outcome of the fitting. This parameter 
can, however, be included in the proposed procedure either after 
the last step by re-building the echelle diagram with a model with 
an adequate Z or by using a reference grid of models calculated 
with representative metallicities. 

5. Preliminary comparison with observations 

Having tested our theoretical approach on artificial pulsators, we 
may now try to evaluate how it can be applied to real obser- 
vations. The difficulty with ground-based observations of PMS 
6 Scuti stars is that these do not provide yet a large set of fre- 
quencies for each star. But space observations have already been 
repo rted where severa l tens of frequencies with very high preci- 
sion ( lMatthewsll2007l) are measured for a single star These will 
be the ideal cases for applying the approach described above. 



Fig. 7. The same as Fig.|6] but for Star!. 

corresponding constant Av lines with the evolutionary tracks are 
inside the assumed uncertainty box (see Figs|4]-(5]l. For the mod- 
els located at these intersections we use the computed frequen- 
cies to construct theoretical echelle diagrams and search for the 
one that best reproduces (i.e. matches as many as possible of) 
the observed frequencies. As a result we obtain M - 2.0 Mq, 
T^s = 8184 K for Starl and M = 2.4 Mq, = 7739 K for 
Starl. The corresponding echelle diagrams are shown in Figs|6]- 

El 

4.4. Results 

As shown in Fig.|6] for Starl all the frequencies are reproduced 
with our best fit model, and we are able to obtain the correct 
spherical harmonic degree / and radial order n, for each fre- 
quency. Our best fit model has also the right mass and an effec- 
tive temperature that differs from the true one by only 200 K. We 
remind that this is the case with a poorly determined luminosity, 
but whose impact on the mass estimation is weak. Based on the 
mode identification provided by this approach, a more detailed 
study of the global characteristics and the internal structure of 
the star could now be attempted. 

Also for Star2 we are able to reproduce the right mass and ef- 
fective temperature and all the 16 frequencies with the right I and 



5.1. The case of V351 Ori 

As a preliminary application we considere d the known multi - 
periodic PMS S Scuti star V351 Ori (see iRipepi et al.ll2003h . 
The observed frequencies for this object, with th e associated 
uncertainties, are reported in Table |4] As noted by 'Ripe pi et al.l 
(|2003|) the first four frequencies are well established, whereas 
/s is slightly less reliable. Moreover and h are quite close to 
/i and fi respectively. As discussed bv lBreger & Bisch^ (|2002|) 
close frequency pairs should be investigated accurately in or- 
der to establish whether they are real or not. In particular the 
quoted authors suggest a method to discriminate which hypoth- 
esis (close pair or single mode with amplitude variation) is cor- 
rect. Unfortunately, current data on V35 1 Ori are insufficient to 
follow this approach. Consequently we will mainly focus on re- 
producing frequencies /i , /2 and . 

The ranges of luminosity and effective temperature as re- 
ported in the literature for V351 Ori are plotted in the HR di- 
agram (dotted box) in Fig. |8] with the PMS evolution tracks. 
Similarly to the plots presented for the two test stars, we also 
report in the same diagram the lines at constant large frequency 
separation. This comparison allows us to obtain a range in mass 
and large separation consistent with the estimated HR diagram 
position of the star, namely 1.8 < M/Mq < 3.0 and 20 fiHz < 
Av„ < 50 yuHz. 

The next step is to identify from the frequency data if a 
stronger constraint can be posed on the large frequency separa- 
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Table 4. Observed frequ encies and associat ed uncertainties for 
V351 Ori, as reported bv lRipepi et"!!] (|2003|) . 





Frequency 


error 

ytiHZ 


/l 


181.6 


0.4 




165.9 


0.4 


h 


147.6 


0.4 




183.8 


0.4 


k 


148.3 


0.4 




4.2 4 3.8 3.6 



Fig. 8. The same as Fig.|4] but for stars V351 Ori (box with dot- 
ted square Une), IP Per (box with dot-dashed line). 



tion. The large separation of V35 1 Ori is difficult to estimate di- 
rectly from the power spectrum. Instead, we use the separations 
between observed frequencies to evaluate for what conditions 
the modes with larger amplitude (expected to be unstable low- 
order radial modes) are consistent with the range of Av found 
from the HR diagram. By using the reference grid we can ex- 
trapolate to the asymptotic regime to find that the frequencies 
are consistent with a Av^ ~ 40 //Hz. 

Consequently, to find the best model that reproduces the ob- 
served frequencies we built the echelle diagram by varying the 
mass and the large separation within the constraints indicated 
above. In doing so we assume that the mode with the largest 
amplitude (/i) is a radial mode. As a result, we find that the 
model in the grid consistent with the observed frequencies has a 
mass M - 2.0 Mq and an effective temperature of T^s - 7539 K. 
As shown in the echelle diagram (see Fig.|9]l, the frequencies /i 
and /3 are associated with radial modes (/ = 0), namely to the 
first and second overtone respectively, while /2 is a nonradial 
mode (rather close to the 1-2 theoretical sequence). As for 
and /s their position in the echelle diagram is between the I - 
and the I - 1 sequences, but as discussed above and in particular 
for /s, the true nature should be investigated carefully on the 
basis of more extensive data. 

There are key assumptions that were made that require fur- 
ther verification. One is the assumption that the highest ampli- 
tude mode (fx) is a radial mode. This is required to properly 





400 










350 


' ' \ 


=0 V 


41/, = 39.5 (jiHz) J 
vo - 25 (mHz) 




300 


1 


^ \ 1-1 


I'o.obs = 30 (^Hz) r 


3, 


250 


















+ 
.is: 


ZOO 








< 


150 










100 










50 









20 40 60 

l'„,l-(Al'a*k + l'o) (a^Hz) 

Fig. 9. The same as figured for V351 Ori. 



adjust the echelle diagram to the observed spectrum. This hy- 
pothesis, or an equivalent calibration for a radial mode, requires 
verification from spectroscopic mode identification. 

Rotation also introduces an additional difficulty as it may 
change the mode identification we obtain for non-radial modes 
and shift the frequencies. In particular, for our best-fitting model 
the rotational splitting should be of the order of 8 mj sin / //Hz 
for an e stimated rotation velocity vsin/ ^ 100 km/sec (see 
iBalona et al. 2002). Such a relatively large splitting removes the 
possibility that the separation between fn and fx and between 
and fs is due to rotation, once the hypothesis of radial modes 
for /i and /3 is abandoned. However, /2 is shifted by about 
10 juHz from the predicted frequency, which may correspond to a 
m mode, pending spectroscopic confirmation. It should also 
be noted that for such a relatively rapid rotation higher-order 
effects must be taken into a ccount (e.g., 'Gough & Thompsonl 
[1990; Dziemb owski & Good e. 1992; Soufi et al., 1998). These 
lead to non-uniform splitting and shifts of modes with m = 0, 
including the radial modes. It is likely that such rotational fre- 
quency perturbations will be a substantial complication in the 
analysis of oscillations of PMS stars, which typically are rapid 
rotators. 

5.2. The case of IP Per 

Another case that we may also consider as a possible prelimi- 
nary application is I P Per, which has be en the target of a recent 
multisite campaign (iRipepi et al.Ll2006h . The observed frequen- 
cies are listed in Table|5]while in Fig.[8]the HR diagram includes 
the location of the star (box with dot-dashed line) and the refer- 
ence grid of evolution tracks. From this plot we can restrict the 
mass and large separation expected for this star to the ranges 
1.6 < M/Mq < 2 and 50 yuHz < ^vt < 75 //Hz. 

Due to the higher number of modes, and considering that 
five of these have a large amplitude, the separation between ob- 
served frequencies indicates that the data are consistent, within 
the expected range, with Av^ ~ 50 //Hz, when extrapolated to 
the asymptotic regime. Following these constraints we have then 
used the echelle diagram to find the best stellar mass and age 
that reproduce as many as possible of the observed frequencies. 
As a result we find a best fit model with M = 1.8 Mq, and 
T^eif = 7773 K. The corresponding echelle diagram is shown in 

According to this plot, frequencies from fx to fn seem to 
align along the sequence for I - modes, whereas /s, and 
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Table 5. Observed fr equencies and assoc iated uncertainties for 
IP Per, as reported bv lRipepi etsi] (l2006h . 





Frequency 
yuHz 


error 

ytiHZ 


/l 


264.9 


1.3 


fl 


400.5 


1.3 
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352.4 


1.3 


/4 


558.2 


1.3 




333.2 


1.3 


h 


277.6 


1.3 
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1.3 
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1.3 
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Fig. 10. The same as figure |6] for IP Per. 

/g are better in agreement with I - 1 theoretical predictions, and 
fg could be an Z = 2 mode. The discrepancy between the the- 
oretical and observed frequencies is probably due to a residual 
uncertainty on stellar mass, intrinsic to the adopted model grid. 
A slightly lower mass would probably allow us to better repro- 
duce the observed frequencies in the echelle diagram. In any case 
the frequency fi is not consistent with a p mode as it lies well 
below the expected frequencies for the range of stellar masses 
of IP Per. This frequency could correspond to a g mode, but we 
have not explored this possibility further here. 

6. Conclusions 

In this work we establish a procedure based on an extended grid 
of models and oscillation frequencies to analyse the properties 
of PMS 6 Scuti stars. The procedure uses a grid of PMS evo- 
lution models to identify in the HR diagram the range of stel- 
lar mass. From the frequencies we estimate the large frequency 
separation which is then used to reduce the uncertainty on stellar 
mass. Finally, from the detailed analysis of the echelle diagram 
for the few possible combinations of mass/age, a fit of as many 
observed frequencies as possible is obtained. This fit provides 
an initial guess on the mode parameters and a precise estimation 
of the stellar mass. The underlying principle of this approach is 
to extract sequentially the information from the frequencies that 
are more robust and directly connected to the global parameters, 
reducing in this way as much as possible the effect on the infer- 
ences of the unknown physics in this type of stars. We also notice 
that the proposed method is quite general, as it can be appUed by 
using models computed with any evolutionary code. 



In order to evaluate the validity of the proposed approach we 
have considered two test cases. In both cases the result was posi- 
tive, and the parameters were recovered with considerable preci- 
sion even when the luminosity of the stars is poorly known. The 
consistency between the true and the inferred stellar parameters 
for both star 1 (computed with STAROX) and star 2 (computed 
with CESAM), seems to indicate that the method does not de- 
pend on the adopted evolutionary code. A comparison with other 
evolutionary codes (e.g. FRANEC, ATON, Siess et al. 2000) will 
be performed when a large number of accurate frequencies, as 
expected from CoRoT observations, will be available for PMS 6 
Scuti stars. 

The method we propose is mainly aimed at studying stars 
whose pulsation data include several tens of frequencies, as it is 
expected to be obtained by CoRoT in the near future. However, 
we have also considered the possibility of applying it to stars 
whose frequency spectrum includes a smaller number of mea- 
sured frequencies. To do so we have reported the preliminary 
applications to two stars that have been targets of multisite cam- 
paigns: V351 Ori and IP Per The uniqueness of the solution 
is strongly dependent on the number of frequencies being used 
when these are less than about ten. However, the results illustrate 
the capability of this approach to constrain the stellar parameters, 
and in particular the stellar mass, even when a small number of 
frequencies is available. 

For a higher number of frequencies the stellar mass and age 
can be constrained and the frequencies used to test the physics 
of the models. The application to a much larger number of fre- 
quencies for PMS pulsators (in particular as expected from space 
observations) and the adoption of finer model grids could consid- 
erably improve our knowledge of the stellar properties and struc- 
ture in this very important phase of stellar evolution. A compar- 
ison with other evolutionary codes (e.g. FRANEC, ATON, Siess 
et al. 2000), as planned in the context of the CoRoT/ESTA col- 
laboration, is also desirable. 
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